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ON THE ACTUAL AND PROBABLE ERRORS OF INTERPO- 
LATED VALUES DERIVED FROM NUMERICAL TABLES BY 
MEANS OF FIRST DIFFERENCES. 



BY E. S. WOODWAED, C. E. 

§ 1. All numerical calculations dependent on tabulated values of loga- 
rithms of numbers, sines, tangents etc., or on natural trigonometric functions 
are subject to certain errors arising from unavoidable inaccuracies in the ta- 
bular values themselves. The inaccuracies in the tabular values result from 
the neglect of figures beyond the tabular number of places, the tabular val- 
ue being given always to the nearest unit in the last place. The actual er- 
ror wliich may exist in the final result of any computation will depend on 
the form of the computation, or the manner in which the possible actual 
errors of the tabular values used are combined. The investigation of this 
class of errors becomes a matter of importance as well as interest whenever 
the highest degree of precision is sought in computations dependent on 
values derived from a given table. Thus, the practical question may arise, 
"What are the limits of accuracy attainable with a 5-place table of loga- 
rithms?" A complete discussion of the various problems that may arise in 
this connection would form an important branch of the theory of errors, and 
might well form the subject matter of a special treatise. In the present 
paper, however, it is proposed to consider only the most elementary of the 
problems, viz., that relative to the error which may exist in an interpolated 
value derived by means of first differences from consecutive tabular values. 
In other words, it is proposed to give some answer to the practical question, 
"What are the possible actual errors and what are the corresponding proba- 
ble errors to which values interpolated from the ordinary 5-place or 7-place 
tables of logarithms are subject?" 

Before taking up this problem, some remarks on the actual and probable 
errors of tabular values are required. 

§ 2. The possible actual errors of tabular values are confined between 
the limits +| and — J of a unit in the last place of a tabular value, and all 
errors between these lemits are equally likely to occur. These facts may be 
taken as self-evident and as a necessary result of the method of constructing 
tables. The probable error, being defined as that error which is as likely 
to be exceeded as not, will be, in such a system of errors, one-half the max- 
imum error, or ±1 of a unit in the last place of a tabular value. 

These statements may be readily verified in any special case. For exam- 
ple, if a large number of actual errors of tabular values of 5-place logari'ms 



—144— 

be derived by means of a 7- or higher place table, it will be found that 

1st, the numbers of positive and negative errors are sensibly equal; 

2nd, the numbers of errors, without regard to sign, lying between equi- 
distant limits are sensibly equal ; and hence 

3rd, one-half the errors will be greater and one-half less than J. 

That the probable error in a system of errors of constant probability be- 
tween given limits is one-half the max. or limiting error may also be proved 
in a more general manner, and some useful principles will result therefrom. 

For this purpose, let s be the general symbol for an error lying between 
the limits ±a. The facility with which e can occur may be denoted by f(s) 
which must be in this case a constant, since by hypothesis all values of e oc- 
cur with the same facility. Hence we may write 

^(s) = c, a constant. 

Supposing £ to vary continuously between the given limits we may make 

(p{e)de = i (is = 1. 

-a ^ —a 

This is the analytical expression of the fact that the sum of all the prob- 
abilities of all possible errors within the given limits is unity. The value 
of the constant c results at once by integration. Thus 

de = 2ao = 1, whence 



^ —a 



Whatever form ^(e) may have, as long as ^(-{-s) = ^( — e), or equal posi- 
tive and negative errors are equally likely to occur, the probable error is the 
limit r in the deiinite integral 

/— r /»0 /•+>■ /•+" /'+<» 

(p{s)ds = i <p{s)ds = I f{e)ds = | ^(e)^ = J t f{e)ds. 

(f{i)dt == J, and ^(e) = --. Hence 

-a ■^C' 

f^^<f{e)de == f^^ 1^ = J, whence r = |. 

§ 3. The errors to which interpolated values are subject may now be 
considered. 

Let v^ and v^ be two consecutive tabular values and v' an interpolated 
value t tenths from v^. Then the value of v' actually computed is 

V' =Vi + <K ~»i) = (1 — i)«i + ^2- (2) 

Denoting the actual errors of Vj, v^ and v' by Sj^, e^ and s' respectively, 

we have 

e' = (1 — i)ei + fej. (3) 
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The meaning of these equations should be clearly understood. Equation 
(2) means precisely what the symbols indicate, and does not imply any neg- 
lect of decimals beyond the tabular number. However, so far as errors 
here considered are concerned, it is practically sufficient to retain one figure 
in addition to the tabular number, retaining for example the 6th figure in 
interpolated values from a 5-place table and the 8th figure in interpolated 
values from a 7-place table. It will be noticed that (vj — '^i) is the tabular 
difference. Equation (3) shows that the error s' of v' cannot exceed ± J, 
since s^ and t^ cannot surpass ± J. An important property of e^ and s^ 
is their mutual independence; i. e., one of them cannot be expressed in 
terms of the other. They are also independent of t and continuous between 
the limits ± |. Therefore e' may have any value between the limits -\- J 
and — ^ for any value of t. 

The species of interpolated value just defined is that which, as will pres- 
ently appear, should be used always if the highest precision with a given 
table is desired. The species of interpolated value more frequently used, 
however, is that formed by abridging the product t{v2 — 'J'l) to the nearest 
unit in the last place of the tabular value. This neglect of decimals is a 
source of error additional to those arising from the errors of the tabular val- 
ues, and may be denoted by Eg. Calling the actual error of the interpolated 
value in this case e", the meaning of Sj^ and Sg being the same as defined 
above, 

e" = (1— Qei+ie^+ss- (4) 

It will be observed that eg in this equat'n differs from s^ and e^'va two 
important respects; viz., 1st, £3 is not independ't of t; 2nd, £3 is disoontin- 
uous for any value of t. Thus, suppose i = J ; then the only possible values 
of Sj are 0, + J and — J. Again, suppose t =^; then the only possible 
values of Sg are 0, + J and — ^. In a similar manner the possible values 
of £3 may be assigned for any value of t The maximum possible value 
of e" varies, therefore, with t A few cases may be written down. 

For i = I or f, the max. of e" = f XJ + ^XJ + | = A- 
" i = 1 or I, " " " £" = fxj + iXi + 4 = 1 . 
" i = Jorf, « " " e" = fXi + IXJ + i = f. 

It will also be noticed that for certain values of £3, s" will be necessarily 
positive or negative whatever values £j and £3 may have. Thus for t = ^, 
if the tabular difference {v^ — tj^) be odd, £3 must be either +-| or — ^, and 
according as £3 is taken as +J or — ^ so will e" be either positive or neg- 
ative, being confined between the limits and -f 1 in the one case and 
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and — 1 in the other. On the other hand, the value of Sg may be such as 
to permit ))oth positive and negative values of e" but with unequal limits 
for the two classes of errors. Thus for i = |-, Sg must be either 0, + i Of 
— J. If £3 = + J, the possible negative values of s" must lie between 
and — ^, while the positive values of e" must lie between and +f • 

Since in a given numerical table the differences must be limited, the pro- 
duct t{v2 — «i) must vanish for < = 0; and hence £3 = for t = 0. 

Besides the two species of interpolated values already defined there may 
be two more, both of little practical importance, but of some theoretical in- 
terest from the fact that their definition gives some light on the real nature 
of the first two species. It is plain that instead of the tabular difference 
{v^ — Vj) there might be used the true differences, such as would be obtained 
practically by inserting in a given table differences derived from a higher 
place table and carried some decimals beyond the last place of values in the 
given table. Would there be any advantage or disadvantage in the use of 
such differences? There may be two cases. Let J be the true difference, 
v-^ the adjacent tabular value and t the interpolating factor. Then, suppos- 
ing no figures omitted, the interpolated value becomes 

^1 + t^. (5) 

The actual error of this value must be the same as the actual error of v-^, 
since there is no error in t^. Denoting for the sake of distinctness the act- 
ual error of «j-f-iJ by e'", we have 

£'" = £1. 

If, however, the product tJ be abridged to the nearest unit of the last 
tabular place, another source of error will be introduced. Calling the actu- 
al error from this source £4 and the actual error of the Interpolated value 
£*", there results 

£- = £, -f £4. (6) 

In this equation e^ may have any value between the limits -f J and — ^ 
for any value of t greater than the reciprocal of twice the greatest difference 
in the table used. 

For values of t greater than this limit £" may have any value between 
the limits -f-l and — 1. For values oft less than the above limit, £4 must 
be less than -f -J, and hence e'^ must be less than ± 1. When t = 0, 
tJ vanishes and hence £4 = 0. 

The characteristics of the four different species of interpolated values and 
the expressions for their actual errors may now be collected in tabular form 
for the purpose of ready comparison. 
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No. of 
species. 



1 

2 
3 
4 



Obtained by 
means of 



Tabular diff. 
a <( 

True differ'e. 
« it 



Expression for actual 
error. 



£"=(l-i)£i+fe2+£3 



Remarks. 



Sj and Ej are indep. of 
< and of each other & 
continu's between ± J. 
Ss is dependent on i & 
discontinu's. e^ is con. 
bet. limits depen't on t. 



§ 4. Thus far the actual errors only of interpolated values have been 
considered. These errors may indeed be ascertained for any special case 
by a suitable computation, but such computation is generally inexpedient 
and of little interest. The probable error of each species of interpolated 
value is of interest and importance, however, since it alFords a measure of 
the precision of such value. The probable error is, as already defined, that 
error which is as likely as not to be exceeded. It applies therefore to the 
aggregate of all possible errors in a system of errors, and may be taken as 
the representative error of the system. It applies in general only to contin- 
uous errors and is itself always a possible error of the system of errors. 

This premised, let r be the general symbol for the probable error in a sys- 
tem of errors of wliich any one is e, and let f{s) express the law of facility 
of these errors. Then, supposing equal positive and negative errors of equal 
facility, or f(-{-e) = <p{ — s), r is defined by the definite integral 

/-r /"O r+r c+i r+i 

(p(£)d£ = I w{£)ds = I (p{e)de = I w{s)ds = | I f(s)de* (7) 

in which ±lia the limiting value of s. The determination of the probable 
error corresponding to s', e", e'" or e'", requires therefore a knowledge of 

It will now be shown what forms ^(e) has for the system of errors rep- 
resented by 

Put (1 — t)si = X and ^2 = y. Then £' = x -j- y, in which x may have 
any value between the limits ±(1 — t)^=a, say, and y may liaveany value 
between the limits ±t^ = b, say. All values of ic and all values of y be- 
tween these limits are equally likely to occur, which facts, as shown in § 2, 
are expressed analytically by 

■*This permits i^{e) to be of such a form that ^(f)de will represent a relative probability. 
If <l>{e)de represent an absolute probability then must 
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<p{x) = ^ and <p{y) = ^. 

The probability of the occurrence of any particular x Is 

<p(x)dx = -r-. 
2a 

Likewise the probability of the occurrence of any particular y is 

Since x and y are independent, the probability of their concurrence is 
f{x)(p{y)dxdy = -^ = -^-^, since e = x^y. 

But since there must be in general an indefinite number of pairs of values 
of a; and y=s — x, which will produce e, (p{£)d£ or the probability of e must 
be equal to da;(?s-=-4a6 taken as many times as there are units in the range 
through which x may vary. In other words, 



<p{B)dt = dej^ 



Aab 

In evaluating this integral, x must not surpass ±a and e-x=ymnst not 
surpass ±6. For any value of e lying between — (a-{-b) and — (a — b), as- 
suming a^b, the limits of the integral are — a and (e + b). This fact is 
rendered plain by a numerical example. Thus, suppose a == dz5 and b = 
±3. Then — (o4- 6) = — 8 and — (a — 6) = ~ 2. Suppose e = — 6 a 
number intermediate to — 8 and — 2. Then the following are the possible 
integer values of x and y which will produce £ = — 6, 

s' X y Limits of x. 

— 6 = — 5--1 — 5 = — a and 

= —4—2 —3 = £'+6 

:= —3—3 

Again, suppose s' is -3, a number between the limits -8 and -2. Then 

e' ^ y Limits of a;. 

— 3 = —5+2 —5 = —a. 

= -4+1 
= —3+0 
= —2—1 
= —1—2 

= 0—3 =: e' + &. 

Therefore, when e lies between — (a+6) and — {a — 6) 
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(8) 



-a ^ab 4:ab ' ^^ ' 4ab 

Similarly, when e lies between — (a — b) and +(a — b), the limiting values 
of a; are e — 6 and e+6. Therefore 



<p{.)de = def-^^^ = ^de, or ^(e) 



(9) 



__ _ -^ 

-6 4a6 4a6"'' "" ^^'^ 4a6' 

Finally, when £ lies between (a — 6) and (a+6) the limits of x are £ — 6 

and +a, and hence 

— s-{-a-{-b 



<p{s) = 



Aab 



(10) 



The "probability curve" represented by equations (8), (9), (10), is made 
up of three straight lines. Equation (8) represents the line joining the p'ts 
whose coordinates are — {a-\- b), and — (a — b), (1 -^-2a). Equation (9) 
represents the line joining the points whose coordinates are -{a-b), (l-^2a) 
and -\-{a — 6), (l-j-2a). Equation (10) represents the line joining the p'ts 
+ (a — b), (1 -i- 2a) H- (a + b), 0. In the special case under consideration, 
a = (1 — ()| and b = t.\. Therefore 

a + b = h 

a-b = i(l-2«), 

L = _L- 

2a ~ l — t' 

The figure shows the "probability curves" cor- 
resp'ing to several values of t. The origin is at 0. 
AOB is the axis of £ and OC the axis of <p{s). 

It may now be observed that eij's (8), (9), (10) 
give 



/ 



4a6 



■ (a+h) 



de+J- 



J ^,„ ^x Aah 



-(<-6)4a6 



ds 



Therefore, by (7) the probable error in the sys- 
tem of errors defined by (8) - (10) is the limit r 
which will satisiy one of the two following con- 
ditions : 

■r £+0+6 ^^ _ p(a+>) _e + a+6 ^^ ^ ^ 
,(a+6) 4a6 J ^^ Aab ^' 

(«) 

n 26 , r+' 26 , , ,a. 



f 



Aab 

[To be continued.] 




